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Abstract 

In this work, we study the radiative leptonic decays of B~ and D~ mesons in the standard 
model and the two-Higgs-doublet-model type II. The results are obtained using the factorization 
procedure, and the contribution of the order 0(Aqcd I^q) is included. The numerical results 
are calculated using the wave-function obtained in relativistic potential model. As a result, the 
decay mode B —)• yrn,- is found to be sensitive to the effect of the charged Higgs boson. Using the 
constraint on the free parameters of the two Higgs doublet model given in previous works, we find 
the contribution of the charged Higgs boson in the decay mode B —)• yrr'r can be as large as about 
13%. 
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I. INTRODUCTION 


Q 


Q. 


A resonance at 126 GeV is discovered in both the ATLAS [l| and CMS experiments 
which is consistent with the Standard Model (SM) Higgs boson j^. This is one of the great 
experimental achievement in recent years. Despite the great success of SM, there are still 
problems implying new physics beyond SM. One of the simplest extensions is the two Higgs 
doublet model (2HDM), with a second Higgs doublet introduced by Ref. jd] (see also js]). 
The Higgs potential and the Yukawa Lagrangian corresponding to the 2HDM is not unique. 
The one with a CP-invariant potential, and with the Yukawa Lagrangian such that one of 


the Higgs doublet couples to the c 
which is so called 2HDM-Type-H 


own sector of fermions while the other couples to the up, 
g. Such a Lagrangian can avoid tree-level flavor changing 
neutral currents (FCNCs) (7|. There are 2 Vacuum Expectation Values (VEV), denoted as 
Ui, and V 2 in 2HDM, and 5 Higgs bosons denoted as h. A, and H^. By studying the 
phenomenons induced by 2HDM, and comparing the results with the experimental data, 
one can detect the signal of the existence of new physics, and hx the free parameters in the 
model. 

The radiative leptonic decay of the heavy pseudoscalar meson with a massive lepton 
provides a good opportunity to study the 2HDM. This process can be mediated both by 
W bosons and bosons, however, the contribution of the charged Higgs bosons is barely 
investigated before. The contribution of 2HDM is usually small. To probe the signal, 
sufficient precision of the numerical results in the SM is required. In this decay mode, the 
strong interaction is involved only in the initial hadronic state. As a result, this decay mode 


B3, 


and results has been calculated 


has been investigated using the factorization ^proach 

up to the order 0{as^^QiCT)/mo) in the SM {q]. Apart from that, with a massive lepton in 
the hnal state, the results are expected to be sensitive to the charged Higgs bosons. 

There are two free parameters in 2HDM associated with this decay mode, tan/? = V 2 / Ui, 
and Mh±- The branching ratios can be expressed as the function of R = fan (3 /Mh± ■ Using 
the constraints on those parameters inl-ll3| , the branching ratios including the contribution 
of the charged Higgs bosons can be obtained. The results are shown in Table IIVI We hnd 
that, the decay mode B —)■ yrz/,- is very sensitive to the charged Higgs bosons, which makes 
it a good decay mode to probe the signal of 2HDM. On the other hand, the decay mode 
D —)• is barely affected by 2HDM, and should be excluded in the search of the charged 
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FIG. 1; tree level amplitudes in SM, the double line represents the heavy quark propagator but 
not the HQET propagator. 

Higgs bosons. 

This paper is organized as follows. Sec. II is a briefly review of the factorization in the SM. 
The contribution of the 2HDM is discussed in Sec. III. The numerical results and analysis 
are contained in Sec. IV. Sec. V is a summary. 

II. FACTORIZATION IN STANDARD MODEL 

The heavy pseudoscalar meson is constituted with a quark and an anti-quark, and one 
of the quarks is a heavy quark. The Feynman diagrams of the radiative leptonic decay of B 
meson at tree level can be shown as Fig. [TJ The contribution of Fig. [T] d is suppressed by a 
factor of 1 , and can be neglected. The amplitudes of Fig. [H a, b and c can be written 


as 



( 1 ) 
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where Pq and pq are the momenta of the anti-quark q and quark Q, p^, pi and Pi, are the 
momenta of photon, lepton and neutrino, e.y denotes the polarization vector of photon, and 
= 7^(1 — 75 ). The idea of factorization is to absorb the infrared (IR) behaviour into the 
wave-function, the matrix element can be written as the convolution of wave-function and 
hard kernel 


14 


F = / dk^k)nUk) 


For simplicity, we denote 


X = rriQ, y = 2/cq • 2 ; = 2p^ - kg, w = 2kQ ■ kq 

Using the wave-function dehned in coordinate space 

^ay{x,y) =< 0\qa{x)[x,y]Qy{y)\f{pq),Q^{pQ) > 


where [x,y] is the Wilson Line which can be written as 15| 


[x, y] = exp 


ips / d^zZfj,A^^{z) 




( 2 ) 


(3) 


(4) 


(5) 


it has been proved that 
can be factorized as 


4 , the matrix element up to the order 0{asAqcD/'mQ) at one-loop 

= J d^kq j d^kq^{kq,kq)Cn{kq,kq,ij,)Tn{kq,kq) 


( 6 ) 


with 


mq^. 


T2 = eq^^Al^P^, T3 = eqP^ " 


T, = eqPl 


2 p'y ‘ kq 

-eh i* 

n = --^P^ + eP\ 


2 kq ■ p^ 

p / g-Pp g -P; 


(7) 


2 kq ■ " 2 p^-pi - “ \p^ ■ Pp p-y-pi, 

except for CiTi, all the other products are contribution of order 0(Aqcd for clarity, 

we dehne Ci = -|- Uj^, with U™ represents order 0((Aqcd/ mg )”*) contribution, the 

coefficients are 
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where yi is defined as: 


X 


log - + -5 + 3 log — + 


yU" x-y 


X , y(2x - 3y) 
(x - y)2 


+ 


log 


X 


y 

( 11 ) 

( 12 ) 


9 ^ / x\ , n X2; , n X , , 

yi = -T^ + 2 Li 2 1-- log^ — + log^ - (13) 

\ yj y y 

In the results given above, the large logarithms at order O ((Aqcd/''^q)°) is already re¬ 
summed. The C 5 T 5 term is the contribution from Ac, and it dose not have 1-loop hard 
scattering kernel, so we have 


6*5 = 1 


(14) 


The numerical result can be obtained by using the wave-function with the one obtained 


5 





























in Ref. 161 




X d-^ikq + k)5-^{kQ - k)6{kqo - \jk^ + 'm?q)5{kQo - yF + tUq) 

with 

Aif = E'’o - '’Q (^.■l'X(-^.t). 'I'® = 4i^mpA5,e-"'-l‘l 

i 

and the result is 

< 7 I qV^^Q \P >= e^ypa£''PpP‘'^{Fv + Fci) 


+ i (e>p • • pp) {Fa + Rci) + Fc2Pp 


with 


Fy = 


1 3 


d^k^{k)- 


{2'kY y/Q J ‘^^yPqoPQoiPqo + mq){pQ0 + uiq) rupE^ 

X {2eqCimQ - C^Pq^Cq + eQ2pqoC3) 


2 ^/PqoPQoiPqQ + rnq){pQ0 + mg) mpE.^ 


X ( 2eqCimQ - C^PgaCq - C 2 - eQ2pgoC3 + Ci 




n -^7 

where Fci and Fc 2 come from the contribution of Ac- Using 


/pP^ = * < Olqr^^QlP >= -t < 0\qril - 75)Q\P > 


we find 


Fci = 


efp p ■ r ,1 ( ^-Pp ^-Pi 
-, Fc 2 = lefpPp 


‘^Pi ■ Pi 


P-1 -Pp P-1- Pi 


and the decay amplitude can be written as 

^*7/VQg 


Asm — 


^/2 


< 7 I qr'‘Q I P > (IP^v) 


(15) 


(16) 


(17) 


(18) 


(19) 


( 20 ) 


( 21 ) 


III. 2HDM CONTRIBUTION 

The Feynman diagrams at tree level of the radiative leptonic decay in 2HDM are shown 
in Fig. [21 The contribution of Fig. IH d is suppressed by Higgs propagator and neglected, 
the decay amplitudes of the others can be written as 


6 



















FIG. 2: Tree level amplitudes of charged Higgs bosons. The double line represents the heavy quark, 
and the dashed line, the charged Higgs. 


A 2HDM _ ieqg'^VQqniQmi tan^ (3 

,2HDM tan2/3 

A =- 0.^2 .^2 -+ 75) 


{p-1 - PqY - ml 




2 Tq • T 7 


(1 + 75)Q(tq) 0 ( 1 - 75 ) 1 ^) 

«(Tq)(^1-75)^) (22) 




2HDM _ -ieg^VqqmQmi tan^ /3 _ 




?(T5)(1 + 75)Q(Tq) 




(1 -75)z/ 


2pz ■ 

where Mh+ is the mass of the charged Higgs boson, and tan /3 is a free parameter in 2HDM. 
The leading contribution at order 0 ((Aqcd/''^q)°) is which will vanish. This can 

be shown immediately, the matrix element can be written as 


< iWqf. 


{Pi -Pqf - ml 


(1 +75)Mq|P > = < '^\Vqf„ 


(P'r - PqY - ml 


” 7'‘(l - 76)^«e|C > 


mq 


(23) 


at the leading order, we can replace pq^/mq by pp^/mp, and the matrix element can be 


factorized as 18 


,ll 7 | 


< l\^qC 


(P7 -Pqf - m\ 


.^M(l _ ry^)uQ\P >= iFAe>^''P'"e^uP'ypPPo 


+ Fc{s3^pp ■ p^ 


p^.e 


7C.7 


■Pp) 

(24) 
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Contracting with pp^, it will vanish. 

In fact, both Fig. [2l a and Fig. [2j b, and the QCD corrections of them will not contribute 
to the decay amplitude because of the symmetry of the wave-function of the pseudoscalar 
mason. Using the similar factorization procedure as in the SM, we hnd that, up to the order 
0(asAQCD/''^Q), the matrix element of Fig. |2J a and Fig. |21 b with 1-loop QCD corrections 
can be factorized as 


Fr, 


2HDM 


(a) = y j I 


d%^{kQ, k,)C™{kQ, kg, kg) 


2HDM/ 


(25) 


the delta-function in the wave-function will replace the transmission momenta kq, kq to the 
on-shell momenta Pq, pq. Using the dehnition of the wave-function in Eqs. flTHll and flTbl) . in 
Dirac representation, we hnd 

F2HDMip) = Yl j d^k^{k)C‘^^^^{x,y,z,w,p)Tr[M ■Tl^^^{pq,pg)] (26) 

n 

with 


y = 2pq-p^, z = 2pg-p^, ,y = 2pq-pg 
^~^)% + ^q)(^ + 70^^75 


M = 


(27) 


\/Pqo{PqO + mq)pqo{pqo + niq) 


and 




2HDM 


= 


7r 7 


2pg ■ p^ 


( 1 + 75)5 Tr 


2HDM 

2 


■yr g 




2HDM 


= e„ 


'yrq 
'2pg-p^ 


(1 + 75 ) 


T|™ = 6 ^( 1 + 75 ) 


2pQ • P7 


i, T5^™ = eQ(l + 75) 


“^Pq ■ Pi 


(28) 


The kinematic indicates that 


p^ = £ = (0,£i,£2,0) 

Pq = (PqO, ki, k2, ks), Pq = {pqo, -ky, -fcs, -^3) 


(29) 


Notice that, are either even functions of ki or ^ 2 , in this case, after the trace, all 

terms are vanished, or the odd function of ki or /c 2 , then, all terms will vanish after the 
integral J d^k. 

However, will survive and will contribute to the decay amplitude. As the same 

case in the SM, also dose not receive contribution from 1-loop hard scattering kernel. 









As a result, we find 




2HDM — 


eGf tan^ {3mimQ 

75 


< 0|,TQ|P > 


+ mi 

2pi ■ 


(1 -75)z/ 


eGf tan^ (3mimQ 




(30) 




T" ■ p^ 


with 


/ 2 HDM — 


1 3 

(27r)3 Ve 


d^km{k) 


‘^^/PqOPQoiPqO + mq){pQQ + TUq) 


[-2mQPqQ - 2pg ■ pq) (31) 


IV. NUMERICAL RESULT AND ANALYSIS OF PARAMETERS OF 2HDM 


The form factors defined in Eq. ffTTD can be calculated using the integral dehned in 
Eq. fflSj) . and we evaluate this integral using [16 1 

mr, = 1.9 GeV, = 5.1 GeV, = 0.08 GeV 

mb = 4.98 GeV, m^ = 1.54 GeV (32) 

Aqcd = 200 MeV, \b = 2.8 GeV-\ Ad = 3.4 GeV^ 

The numerical results of the form factors are inconvenient to use when calculate the decay 
widths. For simplicity, we use some simple forms to ht the numerical results. Inspired by 
the form factors in Ref. 18|], the form factors are htted as 


FAy{E^) — 


AqcD , / Aqcd 




+ B 


Ay 




(33) 


The results are more reliable at the region ^ ^Qcp because we have neglected the higher 
order terms of Aqcd /E^. As we have done in Ref. |9| , we choose the region E^ > 2Aqcd to 
£t the parameters in Eq. fl33l) . The result of B mason in the SM is 


E^{E^) = I 0.25^^^ + 0.39 ( 


E^ 


E^ 


GeV 


-1 


Ff(Ef = 0,; 


A 


QCD 

E^ 


-0.73 


A 


QCD 

E^ 


(34) 


GeV 


-1 


y \ -^7 

For D mason, there is an additional minus sign in E^ and Ey, and the result is 

2 ^ 


E^iE,) = -0.10 


^QCD I 0 76 r 

Ej \ E^ ) 

E^{E^) = fo.39^^^ + 0.04 


GeV 


-1 


(36) 


GeV 


-1 


9 






















0.35 


0.3 - 


0.25 


\ 


> 

Q) 

O 


0.2 - 


0.15 


0.1 - 


0.05 - 


‘. t 






'-k. 


" -K 






-f=- 


0.3 


0.4 


0.5 


0.6 


0.7 


E^(GeV) 


FIG. 3: Fit of the form factors of S —?■ 'yli'i- The points ‘x’ and ‘+’ are the numerical results for 
the form factors Fy and Fa, and the dotted and dashed curves are the fitted results using Eq. p3]l . 
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FIG. 4: Fit of the form factors of D ^ The points ‘x’ and ‘+’ are the numerical results for 
the form factors Fy and Fa, and the dotted and dashed curves are the fitted results using Eq. p3]l . 
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The fitting of the form factors are shown in Fig. [3] and Fig. HI 
On the other hand, Fd and Fc 2 are related to the decay constant, and we use the result 
in Refs. 


16 


18| which is calculated using the same wave-function and the same parameters 


as what we use in this work. The decay constants are 


Jb = 193.57 MeV, fo = 204.98 MeV 


(36) 


The contribution of the 2HDM / 2 hdm can be calculated using Eq. fl3l]) . and the result is 


/ 2 HDM ~ 1-23 GeV , / 2 HDM ~ 0-31 GeV 


D 


With R dehned as i? = tan/? /Mh±, the branch ratios can be written as 


(37) 


Brtot = Brsu x (1 R^niQmia + R^rriQmfb) 


(38) 


with 


a = 


42 


R?mQmi ^g]y[^2HDM + ‘^2HDM'^SM ’ 


h = 




R^rriQUif ^|hdm 


(39) 


Just as the case of the pure leptonic decay 10|, the interference term in the radiative 
leptonic decay is also found to be destructive. We calculate Brsu, a and b separately. Using 
the fitted result of Fa, Fy, the result of fp and / 2 hdm, and using the Gabibbo-Kobayashi- 
Maskawa (GKM) matrix elements 19|, l20|, and mass of the leptons 


Kd = 0.226, Kb = 0.0047, = 1776.82 MeV, m. = 105.658 MeV 


(40) 


we can obtain Brgu, cl and b. There are IR divergences in the radiative leptonic decays when 
the photon is soft or collinear with the emitted lepton. Theoretically this IR divergences can 
be canceled by adding the decay rate of the radiative leptonic decay with the pure leptonic 
decay rate, in which one-loop correction is included 2l|. This is because the radiative 
leptonic decay can not be distinguished from the pure leptonic decay in experiment when 
the photon energy is smaller than the experimental resolution to the photon energy. So 
the decay rate of the radiative leptonic decay depend on the experimental resolution to the 
photon energy which is denoted by AE^. The dependence of the branching ratios of B 
meson in the SM on the resolution are listed in Table HI And for the same reason, a and b 
also depend on AE^, the results of a and 6 of B meson are listed in Table [TTl 
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AE^ 

BrsmiB Br 

Sm(B ^ TVrl) 

AE^ 

BrsuiB pi/^y) Br 

Sm{B TV^i) 

5MeV 

1.35 X 10"® 

1.69 X 10"® 

20MeV 

1.13 X 10"® 

1.27 X 10"® 

lOMeV 

1.24 X 10"® 

1.48 X 10"® 

25MeV 

1.10 X 10"® 

1.20 X 10"® 

15MeV 

1.18 X 10"® 

1.35 X 10"® 

30MeV 

1.07 X 10"® 

1.14 X 10"® 


TABLE I: The branching ratios with different photon resolution AE^ in the SM. 


AE-y 



Jy^B^TU-ry 

AE^ 



JjB^TUry 

5MeV 

-7.47 

14.11 

-2.53 

12.49 

20MeV 

-8.50 

10.01 

-2.47 

10.59 

lOMeV 

-7.98 

12.25 

-2.51 

11.68 

25MeV 

-8.65 

9.19 

-2.46 

10.15 

15MeV 

-8.28 

10.99 

-2.49 

11.09 

30MeV 

-8.77 

8.49 

-2.45 

9.75 


TABLE II: a and b defined in Eq. (j.S9[) with different photon resolution AE^. 


The result oi D ^ tUtI is too small because of the phase space suppression, so we only 
calculate the branching ratios of T* —)■ and the result is listed in Table IIIIl 

Using AE^ = lOMeV 22j, the can be written as 


= 1.24 X 10"® xrH = 1.24 x 10"® x (1 - 4.20/2^ + 3.39i?^) 

Brs^^ru^ = 1.48 X 10“® x th = 1.48 x 10"® x (1 - 22.20R^ + 914.66/?^) 

= 3.64 X 10"® xrH = 3.64 x 10“® x (1 - l.?,2R^ + 0.78i?^) 


(41) 


with Th defined as -Brtot / Br^yi. The numerical results rely on R. In Ref. [10| , two allowed 
regions for R is given as 


0 < Ri < 0.15 GeV-\ 0.22GeV-^ < R 2 < 0.33 GeV^ 


(42) 


In Ref. 


11|. the constraint for Mh+ is Mh+ > 360 GeV at 95% GL, and Ref. 12| shows 


n 


that Mh+ > 315 GeV at 95% GL.. In Ref. [l3(], the tan/? is also constrained. In the case 
that nih = 126 GeV, tan/3 < 5, while in the case that mh < 126 GeV,MH = 126 GeV, 
tan/3 < 30. Gonsidering those constraints, we hnd R < 0.1 GeV, and the R 2 region can be 
excluded. 

The branching ratios including the contribution of charged Higgs bosons in the region 
that 0 < i? < 0.15 GeV“^ are listed in Table HVl And the ratio as a function of R are 
shown in Figs. [5], [6] and [7l 
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FIG. 5; ru as function of R in the B —)• 7^x2/^ decay mode. 



FIG. 6: rH as function of R in the B —>• jTi'r decay mode. 
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AE^ 

BrsuiD 


AE^ 

BrsM{D 


5MeV 

4.42 X 

-7.00 

31.84 

20MeV 

2.87 X 10"® 

-9.47 

25.69 

lOMeV 

3.64 X lO-*^ 

-8.09 

29.46 

25MeV 

2.62 X 10-5 

-9.97 

23.98 

15MeV 

3.18 X 

-8.86 

27.49 

30MeV 

2.42 X 10-5 

-10.38 

22.31 


TABLE III: The branching ratios with different photon resolution AE^ of T> —)• in the SM, 

and the contribution of 2HDM indicated by a and b defined in Eq. (13911 . 



Br in 0 < R < 0.15 GeV ^ 

B 

1.12 X 10-6 < Brtot < 1-24 x 10-6 

B —^ 'Jti't- 

1.28 X 10-6 < Brtot < 1-48 x 10-® 

D 

3.53 X 10-5 < Brtot < 3.64 x 10-5 


TABLE IV: The branching ratios in the region 0 < i? < 0.15 GeV ^ 

We find that, the decay mode B —i- is very sensitive to the contribution of the 

charged Higgs bosons. For R = 0.1 GeV“^, for example, the branching ratios is suppressed 
by more than 13%. 


V. SUMMARY 

In this paper, we calculated the branching ratios of the radiative leptonic decay of the 
heavy pseudoscalar meson with a massive lepton. The contribution of 2HDM-Type-II is 
included. The SM contribution is obtained using the factorization procedure up to the order 
0 (o:sAqcd/ with one-loop correction. The contribution of the charged Higgs boson 
is also obtained by the factorization scheme, however, we hnd that only the diagram with 
the photon emitting from the lepton leg will contribute, which is an order 0(Aqcd/ i^q) 
contribution and dose not receive contributions from 1-loop hard scattering kernel. The 
numerical results of the branching ratios are listed in Table IIVI and dependence of th on R 
is shown in Figs. [5l |6] and [71 

We hnd that, the decay mode B —l- 'JtUt is sensitive to the contribution of the charged 
Higgs in the 2HDM. This decay mode is as sensitive as the pure leptonic decay of B meson, 
which is estimated to be th = {^ — rri^pR?) |^. However, the result of the pure leptonic 
decay is derived at the leading order of 0(( Aqcd/''^q)°), while our result is calculated up to 


14 















FIG. 7: rH as function of R in the D —)■ 7/iz^^ decay mode. 


the order O^asKQCjy/mQ). On the other hand, the branching ratios of this decay mode in 
the SM is about 1.48 x 10“®, which is larger than B —)■ ez/g and B —)■ /iz/^ both are believed 
to be less than 10“® [l9|. The decay mode B -7 y/rz^^ is also important because the charged 


Higgs can suppress the branching ratios as large as about 10%. It is more sensitive than the 
pure leptonic decay modes D —)■ lu, which is estimated to be th = (1 —= (1 —3.6i?^). 
Though the branching ratio of the radiative leptonic decay of the B meson is much smaller, 
the result in the SM is more reliable than the case of the D meson because the heavy quark 
mass rUb is larger then rric- 

We also hnd that, the decay mode D —)■ /rz/^y is rarely affected by the charged Higgs 
boson. On the other hand, the numerical result of this mode is also not sufficiently accurate 
in the SM because the mass of c quark is not so heavy. This decay mode should be excluded 
in the search of the charged Higgs bosons. 
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